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ABSTRACT 



It is usually proposed that hyperaccretion disks surrounding stellar-mass black 
holes at an accretion rate of a fraction of one solar mass per second, which are 
produced during the mergers of double compact stars or the collapses of massive 
stars, are central engines of gamma-ray bursts (GRBs). In some origin/afterglow 
models, however, newborn compact objects are invoked to be neutron stars rather 
than black holes. Thus, hyperaccretion disks around neutron stars seem to exist 
in some GRBs. Such disks may also occur in type-II supernovae. In this paper 
we study the structure of a hyperaccretion disk around a neutron star. Because 
of the effect of a stellar surface, the disk around a neutron star must be different 
from that of a black hole. Clearly, far from the neutron star, the disk may have a 
flow similar to the black hole disk, if their accretion rate and central object mass 
are the same. Near the compact object, the heat energy in the black-hole disk 
may be advected inward to the event horizon, but the heat energy in the neutron 
star disk must be eventually released via neutrino emission because the stellar 
surface prevents any heat energy from being advected inward. Accordingly, an 
energy balance between heating and cooling would be built in an inner region of 
the neutron star disk, which could lead to a self-similar structure of this region. 
We therefore consider a steady-state hyperaccretion disk around a neutron star, 
and as a reasonable approximation, divide the disk into two regions, which are 
called inner and outer disks. The outer disk is similar to that of a black hole and 
the inner disk has a self-similar structure. In order to study physical properties of 
the entire disk clearly, we first adopt a simple model, in which some microphysical 
processes in the disk are simplified, following Popham et al. and Narayan et al. 
Based on these simplifications, we analytically and numerically investigate the 
size of the inner disk, the efficiency of neutrino cooling, and the radial distribu- 
tions of the disk density, temperature and pressure. We see that, compared with 
the black-hole disk, the neutron star disk can cool more efficiently and produce 
a much higher neutrino luminosity. Finally, we consider an elaborate model with 
more physical considerations about the thermodynamics and microphysics in the 
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neutron star disk (as recently developed in studying the neutrino-cooled disk of 
a black hole), and compare this elaborate model with our simple model. We find 
that most of the results from these two models are basically consistent with each 
other. 

Subject headings: accretion: accretion disks — black holes — gamma rays: bursts 
— neutrinos — stars: neutron 



1. Introduction 

Gamma-ray bursts (GRBs) are commonly divided into two classes: short-duration, hard- 
spectrum bursts, and long-duration, soft-spectrum bursts. The observations have provided 
growing evidence that short bursts result from the mergers of compact star binaries and 
long bursts originate from the collapses of massive stars (for recent reviews see Zhang & 
Meszaros 2004; Piran 2004; Meszaros 2006; Nakar 2007). It is usually assumed that both 
a compact-star merger and a massive-star collapse give rise to a central black hole and a 
debris torus around it. The torus has a mass of about 0.01 — 1M and a large angular 
momentum enough to produce a transient accretion disk with a huge accretion rate up to 
~ l.OM s -1 . The accretion timescale is short, e.g., a fraction of one second after the merger 
of two neutron stars, and tens of seconds if a disk forms due to fallback of matter during the 
collapse process. 

The hyperaccretion disk around a black hole is extremely hot and dense. The optical 
depth of the accreting gas is so enormous that radiation is trapped inside the disk and can 
only be advected inward. However, in some cases, this hot and dense disk can be cooled 
via neutrino emission (Narayan et al. 1992). According to the disk structure and different 
cooling mechanisms, flows in the disk fall into three types: advection-dominated accretion 
flows (ADAFs), convection-dominated accretion flows (CDAFs), and neutrino-dominated 
accretion flows (NDAFs). The first two types of flow are radiatively inefficient (Narayan 
et al. 1998, 2000, 2001) and the final type cools the disk efficiently via neutrino emission 
(Popham et al. 1999; Di Matteo et al. 2002). In view of effects of these three flows, 
hyperaccretion disks around black holes have been studied both analytically and numerically 
(e.g., see Popham et al. 1999; Narayan et al. 2001; Kohri & Mineshige 2002; Di Matteo et 
al. 2002; Gu et al. 2006; Chen k Beloborodov 2007; Liu et al. 2007; Janiuk et al. 2007). 

However, newborn neutron stars have been invoked to be central engines of GRBs in 
some origin/afterglow models. First, the discovery of X-ray flares by Swift implies that the 
central engines of some GRBs are in a long-living activity (at least hundreds of seconds) 
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after the bursts (Zhang 2007). This provides a challenge to conventional hyper accretion disk 
models of black holes. Recently, Dai et al. (2006) argued that newborn central compact 
objects in the GRBs could be young neutron stars, at least transiently-existing neutron 
stars, rather than black holes (for alternative models see Perna et al. 2006 and Proga & 
Zhang 2006). These neutron stars may have high angular momentum and their maximum 
mass may be close to or slightly larger than the upper mass limit of nonrotating Tolman- 
Oppenheimer-Volkoff neutron stars. According to this argument, Dai et al. (2006) explained 
X-ray flares of short GRBs as being due to magnetic reconnection-driven events from highly- 
magnetized millisecond pulsars. It is thus reasonable to assume a unified scenario: a prompt 
burst originates from a highly-magnetized, millisecond-period neutron star surrounded by a 
transient hyperaccretion disk, and subsequent X-ray flares are due to a series of magnetic 
activities of the neutron star. 

Second, the shallow decay phase of X-ray afterglows about several hundreds of seconds 
after a sizable fraction of GRBs discovered by Swift has been understood as arising from 
long-lasting energy injection to relativistic forward shocks (Zhang 2007). It was proposed 
before Swift observations that pulsars in the unified scenario mentioned above can provide 
energy injection to a forward shock through magnetic dipole radiation, leading to flattening 
of an afterglow light curve (Dai & Lu 1998a; Zhang & Meszaros 2001; Dai 2004). Recent 
model fitting (Fan & Xu 2006; Yu & Dai 2007) and data analysis (Liang et al. 2007) indeed 
confirm this result. An ultra-relativistic pulsar wind could be dominated by electron/positron 
pairs and its interaction with a postburst fireball gives rise to a reverse shock and forward 
shock (Dai 2004). The high-energy emission due to inverse- Compton scattering in these 
shocks is significant enough to be detectable with the upcoming Gamma-ray Large-Area 
Space Telescope (Yu et al. 2007). 

Third, we note that, in some origin models of GRBs (e.g., Kluzniak & Ruderman 1998; 
Dai k Lu 1998b; Wheeler et al. 2000; Wang et al. 2000; Paczyhski & Haensel 2005), highly- 
magnetized neutron stars or strange quark stars surrounded by hyperaccretion disks resulting 
from fallback of matter could occur during the collapses of massive stars or the mergers of 
two neutron stars. A similar neutron star was recently invoked in numerical simulations 
of Mazzali et al. (2006) and data analysis of Soderberg et al. (2006) to understand the 
properties of supernova SN 2006aj associated with GRB 060218. In addition, hyperaccretion 
disks could also occur in type-II supernovae if fall-back matter has angular momentum. It 
would be expected that such disks play an important role in supernova explosions via neutrino 
emission, similar to some effects reviewed by Bethe (1990). 

Based on these motivations, we here investigate a hyperaccretion disk around a neutron 
star. To our knowledge, this paper is the first to study hyperaccretion disks around neutron 



-4- 



stars related possibly with GRBs. Chevalier (1996) discussed the structure of dense and 
neutrino-cooled disks around neutron stars. He considered neutrino cooling due to electron- 
positron pair annihilation, which is actually much less important than the cooling due to 
electron-positron pair capture in the hyperaccreting case of our interest, since the accretion 
rate assumed in Chevalier (1996) (~ M yr _1 ) is much less than that of our concern (~ 
O.lM s -1 ). In this paper we consider several types of neutrino cooling by using elaborate 
formulae developed in recent years. In addition, we focus on some differences between black 
hole and neutron star hyperaccretion. A main difference is that, the internal energy in an 
accretion flow may be advected inward into the event horizon without any energy release if 
the central object is a black hole, but the internal energy must be eventually released from 
the disk if the central object is a neutron star, because the stellar surface prevents any heat 
energy from being advected inward into the star. Since the accretion rate is always very 
high, the effective cooling mechanism in the disk is still neutrino emission, and as a result, 
the efficiency of neutrino cooling of the entire disk around a neutron star should be higher 
than that of a black hole disk. 

There have been a lot of works to study accretion onto neutron stars in binary systems 
(e.g., Shapiro & Salpeter 1975; Kluzniak & Wilson 1991; Medvedev & Narayan 2001; Frank et 
al. 2002) and supernova explosions (e.g., Chevalier 1989, 1996; Brown & Weingartner 1994; 
Kohri et al. 2005). In the supernova case, spherically symmetric accretion onto neutron 
stars (the so-called Bondi accretion) was investigated in detail. In particular, Kohri et al. 
(2005) tried to use the hyperaccretion disk model with an outflow wind to explain supernova 
explosions. In the binary systems, as the accretion rate is not larger than the Eddington 
accretion rate (~ lO~ 8 M yr -1 ), the physical properties of the disk must be very different 
from those of a hyperaccretion disk discussed here. 

This paper is organized as follows: in §2, we describe a scheme of our study of the 
structure of a quasi-steady disk. We propose that the disk around a neutron star can be 
divided into two regions: inner and outer disks. Table 1 gives the notation and definition 
of some quantities in this paper. In order to give clear physical properties, we first adopt 
a simple model in §3, and give both analytical and numerical results of the disk properties. 
In §4, we study the disk using a state-of-the-art model with lots of elaborate considerations 
about the thermodynamics and microphysics, and compare results from this elaborate model 
with those from the simple model. §5 presents conclusions and discussions. 
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2. Description of Our Study Scheme 

2.1. Motivations of a two-region disk 

We study the quasi-steady structure of an accretion disk around a neutron star with 
a weak outflow. We take an accretion rate to be a parameter. For the accretion rates of 
interest to us, the disk flow may be an ADAF or NDAF. In this paper, we do not consider a 
CDAF. Different from the disk around a black hole, the disk around a neutron star should 
eventually release the gravitational binding energy of accreted matter (which is converted to 
the internal energy of the disk and the rotational kinetic energy) more efficiently. 

The energy equation of the disk is (Frank et al. 2002) 

ds 

Y,v r T— = Q + -Q-, (1) 
dr 

where £ is the surface density of the disk, v r is the radial velocity, T is the temperature, 
s is the entropy per unit mass, r is the radius of a certain position in the disk, and Q + 
and Q~ are the energy input (heating) and the energy loss (cooling) rate in the disk. From 
the point of view of evolution, the structure of a hyperaccretion disk around a neutron star 
should be initially similar to that of the disk around a black hole, because the energy input 
is mainly due to the local viscous dissipation, i.e., Q + = Q^ is - However, since the stellar 
surface prevents the matter and heat energy in the disk from advection inward any more, 
a region near the compact object should be extremely dense and hot as accretion proceeds. 
Besides the local viscous heating, this inner region can also be heated by the energy (Q^ dv ) 
advected from the outer region of the disk. Thus, the heat energy in this region may include 
both the energy generated by itself and the energy advected from the outer region, that is, 
we can write Q + = Qtis + Qtdv- Initially, such a region is so small (i.e., very near the compact 
star surface) that it cannot be cooled efficiently for a huge accretion rate (~ 0.01 — 1M© s _1 ). 
As a result, it has to expand its size until an energy balance between heating and cooling 
is built in this inner region. Such an energy balance can be expressed by Q + = Q~ in the 
inner region of the disk. 

Once this energy balance is built, the disk is in a steady state as long as the accretion 
rate is not significantly changed. For such a steady disk, therefore, the structure of the outer 
region is still similar to that of the disk around a black hole, but the inner region has to be 
hotter and denser than the disk around a black hole, and could have a different structure 
from both its initial structure and the outer region that is not affected by the neutron star 
surface. 

Based on the above consideration, the hyperaccretion disk around a neutron star could 
have two different regions. In order to discuss their structure using a mathematical method 
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clearly, as a reasonable approximation, we here divide the steady accretion disk into an inner 
region and an outer region, called inner and outer disks respectively. The outer disk is similar 
to that of a black hole. The inner disk, depending on its heating and cooling mechanisms 
discussed above, should satisfy the entropy conservation condition Tds/dr oc Q + — Q~ = 0, 
and thus we obtain P oc p 7 , where P and p are the pressure and the density of the disk, and 
7 is the adiabatic index of the disk gas. Also, the radial momentum equation is 

(tf-tf K )r-~(pc 2 s ) = 0, (2) 
p dr 

where Q and Qk are the angular velocity and Keplerian angular velocity of the inner disk, 
and c s = \fP~Jp is the isothermal sound speed. We here neglect the radial velocity term 
v r dv r /dr since we consider the situation v r dv r /dr <C \Q 2 — fiff|r, which can still be satisfied 
when v r <C rQx with Q ~ Qx but \fl — VLk\ > v r /r. Equation ([2]) gives Q oc r~ 3 / 2 and 
c s oc r _1//2 . Moreover, from the continuity equation, we have v r oc (prH)' 1 with the disk's 
half-thickness H = c s /Q oc r. Thus we can derive a self-similar structure in the inner region 
of the hyper accretion disk of a neutron star, 

p OC T - X ^- X \ P OC r-^- 1 ), V r OC r (3-27)/(7-l) ; (3) 

This self-similar structure has been given by Chevalier (1989) and Brown & Weingartner 
(1994) for Bondi accretion under the adiabatic condition and by Medvedev & Narayan (2001) 
and Medvedev (2004) for disk accretion under the entropy conservation condition. In ad- 
dition, if the gas pressure is dominated in the disk, we have 7 = 5/3 so that equation ([3]) 
becomes p oc r~ 3 / 2 , P oc r~ 5 / 2 , and v r oc r -1 / 2 , which have been discussed by Spruit et al. 
(1987) and Narayan & Yi (1994). 

An important problem we next solve is to determine the size of the inner disk. Since 
the total luminosity of neutrinos emitted from the whole disk significantly varies with the 
inner-region size, we can estimate the inner region size by solving an energy balance between 
neutrino cooling and heating in the entire disk. The details will be discussed in §2.3. 

Finally, we focus on two problems in this subsection. First, we have to discuss a physical 
condition of the boundary layer between the outer and inner disks. There are two possible 
boundary conditions. One condition is to assume that a stalled shock exists at the bound- 
ary layer. Under this assumption, the inner disk is a post-shock region, and its pressure, 
temperature and density just behind the shock are much higher than those of the outer disk 
in front of the shock. The other condition is to assume that no shock exists in the disk, and 
that all the physical variables of two sides of this boundary change continuously. We now 
have to discuss which condition is reasonable. 

Let us assume the mass density, pressure, radial velocity, and internal energy density 
of the outer disk along the boundary layer to be pi, Pi, v 1 , and U\, and the corresponding 



-7- 



physical variables of the inner disk to be p 2 , P 2 , v 2 , and u 2 at the same radius. Thus the 
conservation equations are 

P1V1 = p 2 v 2 

Pi + ptvj = P 2 + p 2 v\ (4) 
(ui + Pi + Pl v 2 1 /2)/p 1 = {u 2 + P 2 + p 2 v 2 2 /2)/p 2 . 

From the Rankine-Hugoniot relations, we know that if P 2 3> Pi, the densities of two sides 
of the boundary layer are discontinuous, which means that a strong shock exists between 
the inner and outer disks. On the other hand, if Pi ~ P 2 , we can obtain p 1 ~ p 2 , which 
means that only a very weak shock forms at this boundary layer, or we can say that no shock 
exists. Therefore, we compare Pi and p\v\ of the outer disk. If Pi <C pivf or c s < ui, we can 
assume that a stalled shock exists at the boundary layer, i.e., the first boundary condition is 
correct. If Pi ^> p\v\ or c s ^> Vi, otherwise, we consider P x ~ P 2 , and thus no shock exists. 
In §3 and §4, we will use this method to discuss which boundary condition is reasonable. 

Second, what we want to point out is that the effect of the magnetic field of the central 
neutron star is not considered in this paper. We estimate the order of magnitude of the Alfven 
radius by using the following expression (Frank et al 2002), ta — 2.07 x 10 4 m d 2 ^ 7 m~ 1 / 7 p^ 7 , 
where m d = M/O.OlM s _1 is the accretion rate, m = M/1AM Q is the mass of the neu- 
tron star, and p 30 is the magnetic moment of the neutron star in units of 10 30 Gcm 3 . 
Let r* be the neutron star radius. If the stellar surface magnetic field B s < P SjCr = 
2.80 x lO^mJ^m 1 / 4 ^ 574 G or B s < P SjCr = 0.89 x 10 15 mj 2 m 1 / 4 G when r* = 10 6 cm, 
then the Alfven radius is smaller than tv The critical value B scr increases with in- 
creasing the accretion rate. This implies that the stellar surface magnetic field affects the 
structure of the disk significantly if B s > P SiCr ~ 10 15 — 10 16 G for typical accretion rates. 
Therefore, we assume a neutron star with surface magnetic field weaker than B S)Cr in this 
paper. This assumption is consistent with some GRB-origin models such as Kluzniak & 
Ruderman (1998), Dai & Lu (1998b), Wang et al. (2000), Paczyhski & Haensel (2005), and 
Dai et al. (2006), because these models require a neutron star or strange quark star with 
surface magnetic field much weaker than P S)Cr - 

2.2. Structure of the outer disk 

Here we discuss the structure of the outer disk based on the Newtonian dynamics and 
the standard a-viscosity disk model for simplicity. The structure of the hyperaccretion disk 
around a stellar-mass black hole has been discussed in many previous works. The method 
and equations we use here are similar to those in the previous works since the outer disk is 
very similar to the disk around a black hole. 
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We approximately consider the angular velocity of the outer disk to be the Keplerian 
value £Ik — \/ GM/r 3 . The velocity should be modified in a relativistic model of 
accretion disks (Popham et al. 1999; Chen & Beloborodov 2007), but we do not consider 
it in this paper. We can write four equations to describe the outer disk, i.e., the continuity 
equation, the energy equation, the angular momentum equation and the equation of state. 

The continuity equation is 

M = A7crpv r H = 27rr£t> r , (5) 

where the notations of the physical quantities have been introduced in §2.1. 

In the outer disk, the heat energy could be advected inward and we take Q~ dv = 
( 1/2) TTv r ds / dr to be the quantity of the energy advection rate, where the factor 1/2 is 
added because we only study the vertically-integrated disk over a half-thickness H. Thus 
the energy-conservation equation ([ID is rewritten as 

Q + = Q;, d + Q: dv + Qu- (6) 

The quantity Q + in equation ([6]) is the viscous heat energy generation rate per unit 
surface area. According to the standard viscosity disk model, we have 

Q+ - ^/ (7) 

where / = 1 - (r*/r) 1/2 (Frank et al. 2002). 

The quantity Q~ &A in equation ([S]) is the photon cooling rate per unit surface area of the 
disk. Since the disk is extremely dense and hot, the optical depth of photons is always very 
large and thus we can take Q~ d = as a good approximation. 



The entropy per unit mass of the disk is (similar to Kohri & Mineshige 2002) 



■5gas 



S gas , S vad ^ (5k B k B . \{2Tik B Tfl 2 ]\ 2 



2 aT 



y, 



P 



where the summation runs over nucleons and electrons, k B is the Boltzmann constant, h is 
the Planck constant, a is the radiation constant, Sq is the integration constant of the gas 
entropy, and the term 2g#aT 3 /3p is the entropy density of the radiation with g* = 2 for 
photons and g* = 11/2 for a plasma of photons and relativistic e + e~ pairs. We assume that 
electrons and nucleons have the same temperature. Then we use equation ([8]) to calculate 
ds/dr and approximately take dT/dr w T/r and dp/ dr w p/r to obtain the energy advection 
rate, 

R , . 4 aT 3 " 

{l+Y e ) + -g m — , (9) 



Qadv = V r T— 



2 y c/ 3" p _ 
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where the gas constant R = 8.315 x 10 7 ergs mole -1 K -1 and Y e is the ratio of electron to 
nucleon number density. The first term in the right-hand bracket of equation (Q comes from 
the contribution of the gas entropy and the second term from the contribution of radiation. 

The quantity Q~ in equation (JS)) is the neutrino cooling rate per unit area. The expres- 
sion of Q~ will be discussed in detail in §3 and §4. 

In this paper we ignore the cooling term of photodisintegration Qp hotodis , and approx- 
imately take the free nucleon fraction X nuc rs 1. For the disks formed by the collapses of 
massive stars, the photodisintegration process that breaks down a-particles into neutrons 
and protons is important in a disk region at very large radius r. However, the effect of 
photodisintegration becomes less significant for a region at small radius, where contains less 
a-particlea3. On the other hand, for the disks formed by the mergers of double compact 
stars, there will be rare a-particles in the entire disk. As a result, we reasonably take all the 
nucleons to be free (X nuc « 1) and neglect the photodisintegration process, since we mainly 
focus on small disks or small regions of the disks (as we will mention in §3.4 with the outer 
boundary r out to be 150 km.) 

Furthermore, the angular momentum 
written as 

uE 



conservation and the equation of sate can be 



£/. do) 

37T 



P = P e + P nuc + P Tad + P u , (11) 

where v = ac s H in equation (ITU1) is the kinematic viscosity and a is the classical viscosity 
parameter, and P e , P nuc , -P ra d and P v in equation (fTTj) are the pressures of electrons, nucleons, 
radiation and neutrinos. In §3 we will consider the pressure of electrons in extreme cases 
and in §4 we will calculate the e pressure using the exact Fermi-Dirac distribution function 
and the condition of /^-equilibrium. 

Equations (JSJ), (jSD, (fTO]) and (TITi) are the basic equations to solve the structure of the 
outer disk, which is important for us to study the inner disk. 



1 Kohri et al. (2005), Chen & Beloborodov (2007), and Liu et al. (2007) discussed the value of nucleon 
fraction X nuc and the effect of photodisintegration as a function of radius for particular parameters such 
as the accretion rate and the viscosity parameter a. The former two papers show that X nuc rs 1 and 
Qphotodis ~ for r < 10 2 r 9 . Although the value of X nuc from Liu et al. (2007) is somewhat different from 
the previous works, the ratio of Qph otoc j is /Q + in their work also drops dramatically for r < 10 2 r g . Therefore, 
it is convenient for us to neglect the photodisintegration process for r < 1 2 r g or r < 400 km for the central 
star mass M = 1.4M Q . 
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2.3. Self-similar structure of the inner disk 



In §2.1 we introduced a self-similar structure of the inner disk, and described the method 
to determine the size of the inner disk. Now we will establish the energy conservation 
equation in the inner disk. We assume that f is the radius of the boundary layer between 
the inner and outer disks, and that p, P, and v r are the density, pressure and radial velocity of 
the inner disk just at the boundary layer respectively. From equation ([3]), thus, the variables 
in the inner disk at any given radius r can be written by 

p = p(r/r) 1/{ ' y - 1 \ P = P^/ry'^.Vr = v r (f /r)^ 1 "^'^ . (12) 



We take the outer radius of the accretion disk to be r out . The total energy per unit time 
that the outer disk advects into the inner disk is (Frank et al. 2002) 



^ v . (1 _ A) ™ i 

4 r 



1 - 



1/2' 



r out 



1 - 



r out 



1/2' 



(13) 



where f v is the average neutrino cooling efficiency of the outer disk, 

f!° ut Q-2nrdr 



fv 



jT out Q+2-nrdr 



(14) 



If the outer disk flow is mainly an ADAF, Q~ <C Q + , then /„ ~ 0; if the outer disk flow 
is mainly an NDAF, Q~ ^> Q~ ad and Q~ 3> Q~ dv , we have /„ ~ 1. The heat energy in the 
inner disk should be released more efficiently than the outer disk, we can still approximately 
take Q ~ in the inner disk, and the maximum power that the inner disk can release is 



3GMM 




r out 



1 



2 

3 V 7" out 

f 



1/2 



Tout 



1 



r ont 



1/2 



(15) 



where we have integrated vertically over the half-thickness. The first term in the right-hand 
of this equation is the total heat energy per unit time of the entire disk, and the second term 
is the power taken away through neutrino cooling in the outer disk. Considering r out ^> r*, 
we have 

3GMM { \_ ^ U 
3r* f 



T ~ 

J-'v.max 



4 



_ 2 /r*N 1 - 
3 If 



(16) 



The maximum energy release rate of the inner disk is GMM / Ar* if the outer disk flow is 
mainly an ADAF and f v ~ 0. This value is just one half of the gravitational binding energy 
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and satisfies the Virial theorem. If the outer disk flow is an NDAF, then the energy release 
of the inner disk mainly results from the heat energy generated by its own. 

Following the above consideration, most of the energy generated in the disk around a 
neutron star is still released from the disk, so we have an energy- conservation equation, 

where e is a parameter that measures the efficiency of the energy release. If the central 
compact object is a black hole, we have e ~ and the inner disk cannot exist. If the central 
compact object is a neutron star, we can take e ~ 1 and thus we are able to use equation 
( I17p to determine the size of the inner disk. 



3. A Simple Model of the Disk 

In §2, we gave the equations of describing the structure of a hyper accretion disk. How- 
ever, additional equations about microphysics in the disk are needed. In order to see the 
physical properties of the entire disk clearly, we first adopt a simple model for an analytical 
purpose. Comparing with §4, we here adopt a relatively simple treatment with the disk 
microphysics similar to Popham et al. (1999) and Narayan et al. (2001), and discuss some 
important physical properties, and then we compare analytical results with numerical ones 
which are also based on the simple model. 

If the disk is optically thin to its own neutrino emission, the neutrino cooling rate can be 
written as a summation of four terms including the electron-positron pair capture rate, the 
electron-positron pair annihilation rate, the nucleon bremsstrahlung rate and the plasmon 
decay rate, that is, Q~ = (q eN + q e+e - + g brems + g p i asm on)# (Kohri & Mineshige 2002). 
We take two major contributions of these four terms and use the approximative formulae: 
q e + e - = 5 x lO^T^ergscm" 3 s" 1 , and q e ^ = 9 x lO^pT^ergscm^ 3 s _1 . Thus equation (jSJ) 
can be rewritten as 



3GMM j _ MT 
87rr 3 Anr 2 



R„ 22 aT 3 

— (1 + Y e ) H 

2 v ; 3 p 



(5 x lO 33 !^ + 9 x 10 23 pTn)^-. (18) 



If neutrinos are trapped in the disk, we use the blackbody limit for the neutrino emission 
luminosity: Q~ ~ (|<7bT 4 )/t, where r is the neutrino optical depth. We approximately 
estimate the neutrino optical depth as (q e + e - + q e ^)H/(4: x |(XbT 4 ). 

Moreover, we take the total pressure in the disk to be a summation of three terms P e , 
P nuc and P ra d, and neglect the pressure of neutrinos: P e = n e kT+Ki (pY e )^ 3 , P nuc = n nuc kT, 
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and P ra d = llaT 4 /12, where K\ (pY e ) 4 ' 3 is the relativistic degeneracy pressure of electrons, 
n e and n nnc are the number densities of electrons and nucleons respectively. Here we also 
neglect the non-relativistic degeneracy pressure of nucleons. We thus obtain 

P = P e + P nuc + P rad = p (1 + Y e ) RT + K x ( P Y e f 3 + ^aT 4 , (19) 

where K x = ^ {j^)^ = 1-24 x 10 15 cgs. 

Equations (fTUj) . (fl8|) and (fl9|) can be solved for three unknowns (density, temperature 
and pressure) of the steady outer disk as functions of radius r for four given parameters a, 
Y e , M and M in the simple model. Once the density, temperature and pressure profiles are 
determined, we can present the structure of the outer disk and further establish the size and 
the structure of the inner disk. 



3.1. The outer disk 

We analytically solve equations ( TTUl) . (fl8l) and (fT9l) in this subsection. Our method is 
similar to that of Narayan et al. (2001). However, what is different from their work is that 
we use the same equations to obtain both ADAF and NDAF solutions in different conditions. 
Besides, we find that the factor f — 1 — (r^/r) 1 / 2 cannot be omitted because it plays an 
important role in determining the disk structure. In this subsection, for convenience, we 
expand the solution range to the entire disk (i.e., r* < r < r out ) rather than just consider it 
in the outer region. The size of the inner disk, which depends on the structure of the outer 
disk, will be solved in §3.2. 

First we show a general picture. If the accretion rate is not very high, most of the energy 
generated in the disk is advected inward and we call the disk as an advection-dominated disk. 
As the accretion rate increases, the density and temperature of the disk also increase and the 
neutrino cooling in some region of the disk becomes the dominant cooling mechanism. Thus 
we say that this region becomes neutrino-dominated. When the accretion rate is sufficiently 
large, the disk may totally become neutrino-dominated. Besides the accretion rate, there are 
some other factors such as the mass of the central neutron star, M, and the electron-nucleon 
ratio, Y e , are able to influence the disk structure. 

We take the mass of neutron star M = \AmM & , the accretion rate M = x 
O.OIMqS- 1 , a = 0.1a_i, r = 10 6 r 6 cm, p = 10 n pn gem" 3 , T = 10 n T n K, and P = 
10 29 P29ergscm _3 . In the case that the disk flow is an ADAF with the radiation pressure to 
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be dominated, we find that the density and temperature in the disk are 



T u = 0.832m 1 /8my 4 / 1/8 a:i / V 6 - 5/8 . 



(20) 




(21) 



We have to check the validity of the assumption made in deriving the above solution, i.e., 
we need the relations Q~ dv > Q„ and P ra d > P gas , Prad > Picg to be satisfied, where P ra d and 
Pdog are the gas and degeneracy pressure in the disk. Then we get the relations that 



which means that a larger radius allows a larger upper limit of the accretion rate for the 
radiation-pressure-dominated region. On the other hand, if the parameters m, a and in 
some region of the disk do not satisfy inequations (|22|) or (|25|) . the other types of pressure 
can exceed the radiation pressure but the region can still be advection-dominated. For 
an analytical purpose, we want to discuss the range of different types of pressure in two 
extreme cases where Y e ~ 1 or Y e <C 1. From inequation (I23p . we can see that, since the 
minimum value of r^f' 7 ^ 3 is 19.9, when rhd > 0.453m 7//2 a_i for Y e ~ 1, or rrid > 7.25m 7 / 2 a_i 
for Y e <C 1 , the gas pressure takes over the radiation pressure in the disk and the entire 
disk becomes gas-pressure-dominated. On the other hand, the degeneracy pressure is larger 
than the radiation pressure at a very large radius if the electron fraction Y e is not very 
small. However, we do not consider this situation for an ADAF region of the disk, because 
the degeneracy pressure, even if larger than the radiation pressure, cannot exceed the gas 
pressure. 

When the gas pressure is dominated and the outer disk flow is still an ADAF, we obtain 
the density, temperature and pressure of the disk as 



r 6 / 1/5 > 2.28m~ 3/5 m d /5 aZ 2 v 

r 6 /- 7 / 3 < 74.6 (1 + Y e y 8/3 m 7 / 3 m- 2/3 a 2 J?, 
r 6 f- J / 3 < 17Am 7 / 3 m~ 2/3 a 2 J 3 Y- s / 3 . 



(22) 



(23) 
(24) 




(25) 



r, rr fi /l , TT\-12/11 5/11 • 8/11 /5 /H -8/11 -21/ 

pn = 0.556 (1 + Y e ) ' m 0/ii m d ' j 0/L1 a_ 1 r 6 
Tn = 1.29 (1 + Y e y 3/U m^mf V 4/11 a"f V /U 
P 29 = 5.98 (1 + F e )- 4/n m^W d 0/u f^aZl 0/11 r e 29/ 



(26) 
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Similarly, we discuss the validity of the solution fl26|) . The assumption of a gas-pressure- 
dominated ADAF disk (Q^d v > Q~) requires 

r 47/2 2/ - 20 /ll > 12g (1 + ^-26/1! m 29/2 2rh 10/ll a = 2 l/ll_ (2?) 

-Pgas > Pad can be written by 

r 6 /- 7 / 3 > 74.6 (1 + r e )" 8/3 m 7 ^m d 2/ \ 2 _!l (28) 

-Pgas > -Pdeg leads to 

r 6 r ? / 3 < 8.07 x 10 3 (1 + Y e ) 12 Y^m'^m^a 2 ^, (29) 

which is satisfied for a large parameter space. 

If Y e ~ 1 and the parameters m and a make 0.453m 7//2 a_i < rrid < 1.73m~ 29//20 a 21 / 10 
valid, the entire disk becomes an advection-dominated disk with the gas pressure to be 
dominated. However, if Y e -C 1, such a disk cannot exist, since inequation (1271) cannot 
be always satisfied in the entire disk and some region of the disk would become neutrino- 
dominated. Also, when the mass accretion rate rhd becomes higher, most region of the disk 
becomes neutrino-dominated. 

Furthermore, in the region where neutrino cooling is efficient and the gas pressure dom- 
inates over the degeneracy pressure, we have another particular solution 

Pll = 2.38 (1 + Y e y 9/5 m^m d faZl 3/10 r 6 51/2 ° 

T n = 0.490 (1 + Y e ) 1/b m^a^r- 3/10 (30) 
P 29 = 9.71 (1 + Y e r 3/5 m^m d f a Zl 1/10 r G 57/2 °. 

In this case, the temperature is independent of the accretion rate m^. Finally we check 
the gas pressure-dominated assumption. Using equation (I3"UI) and assuming P gas > Pd eg , we 
obtain 

r 6 r 2(V33 > 3.18r e 80 / 33 (1 + F e )- 36/n m^mf 3 Vf /33 . (31) 

Inequation (|3T|) is always satisfied if Y e <C 1, and thus we can say that the gas pressure- 
dominated assumption is valid. However, if Y e ~ 1, a part of the disk becomes degeneracy 
pressure-dominated if > 64.5ai 9 / 10 m -11 / 20 . m particular, in the case of Y e ~ 1 and large 
accretion rate rhd, a set of solutions on the part of the disk are 

pn = 1.26F e - 4/3 m 2 / 3 mf f^a'Trf 

T u = 0.526F e 4/27 m 13 / 108 my 27 / 1/27 « 7/ 1 5 V 6 - 13/36 (32) 
P 29 = 7.85F e - 4/9 m 8 / 9 m 8/9 / 8 / 9 a- 8/9 r 6 - 8/3 . 
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which describe an NDAF with the degeneracy pressure to be dominated. However, we should 
remember that in deriving the above solutions we have not considered the constraint of r > f. 

Here we make a summary of § 3.1. We used an analytical method to solve the den- 
sity, pressure and temperature of the outer disk based on a simple model discussed at the 
beginning of §3. The accretion flow may be ADAF or NDAF with the radiation, gas or 
degeneracy pressure to be dominated. We fix radius r and show several possible cases in 
the disk with different accretion rate in Table 2. Moreover, for different electron fraction 
Y e and fixing m = 1 and a„i = 1, we calculate the upper limit of m^. If Y e ~ 1, the 
advection-dominated region in the disk can be radiation or gas pressure-dominated, and the 
neutrino-cooled region can be gas or degeneracy pressure-dominated. However, if Y e <C 1, 
the gas pressure-dominated region in the ADAF case is very small, and the degeneracy 
pressure-dominated region cannot exist. In §3.4 we will obtain similar results by using a 
numerical method. 

The solutions given in this subsection can also be used to discuss the properties of the 
disk around a black hole. Our analytical solutions of the outer disk are similar to those of 
Narayan et al. (2001) and Di Matteo et al. (2002), who took Y e to be a parameter. Narayan et 
al. (2001) found that advection-dominated disks can be radiation or gas pressure-dominated, 
and neutrino-dominated disks can be gas or degeneracy pressure-dominated instead. This 
is consistent with our conclusion for Y e ~ 1. However, these authors did not consider the 
factor / = 1 — which is an important factor because a small r, as we have mentioned 

above, can dramatically change the parameter space of the outer disk. Di Matteo et al. 
(2002) discussed different pressure components (their Fig. 2), which is also consistent with 
our conclusion. Chen & Beloborodov (2007) calculated the value of Y e and showed that 
Ye <C 1 when r is small. According to the above discussion, therefore, the degeneracy- 
pressure-dominated region in the NDAF disk cannot exist. This is consistent with Chen 
& Beloborodov (2007) that the pressure in a neutrino-cooled disk is dominated by baryons 
(gas). 

However, our analytical results are partly different with Kohri et al. (2005) and Chen & 
Beloborodov (2007) which showed that the electron pressure is dominant in some advection- 
dominated regions of the disk. This difference is mainly because that we take P ra d = 
llaT 4 /12 in our analytical model, which includes the contribution of relativistic electron- 
positron pairs. However, Kohri et al. (2005) and Chen & Beloborodov (2007) took P ra( j = 
aT 4 /3 and calculated the pressure of e + e~ pairs in the term of electron pressure P e . As 
a result, the radiation pressure we consider in this subsection is actually the pressure of a 
photon and e + e~-pair plasma. 

In the following subsection, we will establish the structure of the inner disk depending 
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on the outer disk solutions here, and use the value of r to further constrain the solutions 
that we have obtained. 



3.2. The inner disk 

3.2.1. Boundary layer between the inner disk and outer disk 

We use the method discussed in §2.1 and compare the radial velocity with the local 
speed of sound of the outer disk using the results given in section §3.1. 

The radial velocity of the outer disk at radius r is 

M M MVGM 

Vl ~ 2vrr£ ~ Aixr Pl H ~ Anr 5 / 2 y/P^' ^ ' 



Hence, we have 



* - 0.04654^. (34) 



c 2 s Pi l6n 2 r 5 P? ' r|P x 2 



2<) 



In the case where the outer disk is ADAF and the radiation pressure is dominated, using 
solutions (J2"0j) and (j2"T|) . we have 

4 ~ 4.22 x 10-Vi/- 1 . (35) 

Therefore, we see V\ <C c s for typical values of the parameters. 

In the case of ADAF with the gas pressure to be dominated, using solution (1261) . we 

have 

% ~ 1.302 x 10- 3 (1 + Y e f n m- 7 / u mf 1 a 20 / 11 r 18/1 V 6 3/11 . (36) 

C s 

For NDAF, using expression (1301 to compare the radial velocity with the speed of sound, 
we have 

4 ~ 4.94 x lO" 4 (1 + Y e f 5 m-^alff-Vj 10 . (37) 

Note that this ratio is independent of the accretion rate. From equations (1361) and (137)1 . we 
still find that v\ c s is always satisfied except for a region very near the stellar surface, 
where the factor / is very small. This region, however, is so small that it belongs to the inner 
disk where we have to use the self-similar structure, which we will discuss later in details. 

Therefore, for the hyperaccretion disk discussed in this paper, as the disk is extremely 
hot and dense, the radial velocity is always subsonic. So there is no stalled shock between 



the inner and outer disks. Thus, all physical variables between two sides of the boundary 
layer between two regions of the disk change continuously. Besides, the rotational velocity 
is assumed to be the Keplerian value and has no jump at the boundary layer. 



3.2.2. Solution of the inner disk 

We now study the inner disk analytically based on the results given in §3.1. The main 
problem that we should solve in this subsection is to determine the size of the inner disk 
for a range of parameters and to describe the structure of the inner disk. In the case where 
the radiation pressure is dominated, by using the self-similar structure (|12|) . we obtain the 
temperature of the inner disk, 

(P \ !/ 4 2-7 - 7 

= LOlm^myV^ali 76 ^ ^ • (38) 

where f and p are the radius and density of the boundary layer between the inner and outer 
disk, and / = 1 — (7%,/f) 1 / 2 . 

Using the self-similar condition and the above expression of Tn, we find the total neu- 
trino cooling rate, 

f r Q-2irrdr = 10 51 x f \2^m^mfaZTf^rP^rf^Hr & . (39) 

The outer disk flow is mainly an ADAF, using the solution of a radiation pressure- 
dominated ADAF (i.e., solutions [201 and I2T1) . we have the total neutrino cooling rate, 

(nj — 1 \ , , 25-157 / 37-8 37-8 \ 

m^mff^aZTr^- 1 " ( r 2 ^ ~ 1 . (40) 

From the energy-conservation equation (fT7|) . the position of the boundary layer satisfies the 
following equation, 

5(5-37) / fZ \ 3 / 4 / 37-8 37-8 \ /O _ Q 0/ \ 

(l - ^f) (rJP - ff^ J = 0.0597 {^1) m^m d 3/2 a^r-l (41) 

where we take e ~ 1 and f v ~ in equation (IT7|) . The left-hand term of equation (T4T1) 
increases with increasing f 6 , so f 6 decreases if rfid increases in the right-hand term of this 
equation. In other words, the size of the inner disk decreases as the accretion rate increases. 
From equation fj4TT) . we can also see that its solution, f^, is independent of Y e , and increases 
with the mass of the central star. In addition, since the gas pressure has its own contribution 
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to the disk, the actual adiabatic index 7 is larger than 4/3, which can makes the solution of 
equation (j4Tj) larger. For an analytical purpose, we assume several different sets of parameters 
to solve equation (j4"Tl) . Table 3 clearly shows that, in the radiation-pressure-dominated disk 
with an advection-dominated outer region, as the accretion rate increases, the value of f 
decreases, and that f increases with increasing 7 or decreasing m. 

In the case where the gas pressure is dominated and the outer disk flow is still an ADAF, 
we obtain the temperature of the inner disk, 



r u = (i + y e - 



1 B 



29 



8.315 p n 

Similarly, from equations ffTTj) and ff26l) . we have 



(1 + Y e )~ l 0.874m 2 / 3 p- 1 2/3 mf / 2/3 alf 3 f 6 ~ V- (42) 



Q-2nrdr = 10 51 x / 45.0p n 10/3 (1 + Y e y 6 m^mf'f^^'^K'^^dr^ (43) 



The neutrino luminosity of the inner disk reads 



37-2 

1 1 3 \ / 2-37 2-37 

w 51/22 • 21/11 ?3l -21/11 -(—-22) 



r 7-1 — r 7-1 
1 *,6 '6 



(44) 



The energy-conservation equation of the inner disk is 

-31/11 



2-37 
, 7-1 
*,6 



2-37 
r 7-1 



382(7 - i; 
37-2 



(45) 



From equation ([431) . we see that the size of the inner disk (f) also decreases with in- 
creasing the accretion rate m^. Table 4 gives solutions of equation fj4"5l) with different sets of 
parameters. We can see that r also decrease with increasing 7, m or decreasing Y e . 

We above study the case where the outer disk is advection-dominated, and find that the 
size of the inner disk always increase with the accretion rate. In the case where the outer 
disk is mainly neutrino-dominated, using expression (130]) and equation fTl7|) . we obtain an 
energy-conservation equation in the inner disk, 



2.77 



' l-l 
.37-2 



27-1 _ / 2-37 

■ ; 1 /['■: 



6 



*.6 



2-37 
r 7_1 



0.924 
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1/2 



(46) 
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This equation shows us that re in NDAF is independent of Y e , m, rhd, and a, but only 
dependent on 7 and f u . The size of the inner disk is constant no matter how much the 
components, the accretion rate of the disk, and the mass of the central neutron star are. 
If the outer disk is mainly an NDAF, we have f v ~ 1. We choose several different sets of 
parameters to obtain the solution of equation (I4"6j) (see Table 5). As f v increases, the value 
of f decreases slightly. Here we also consider an intermediate case of 7 between 5/3 and 4/3, 
the decline of 7 makes the size of the inner disk decrease. 

However, in the above discussion about an NDAF, we have not considered the effect of 
neutrino opacity but simply assumed that neutrinos escape freely. Actually, if the accretion 
rate is sufficiently large and the disk flow is mainly an NDAF, and the disk's region near 
the neutron star surface can be optically thick to neutrino emission. With increasing the 
accretion rate, the area of this optically thick region increases. We now estimate the effect of 
neutrino opacity on the structure of the inner disk. Let the region of r* < r < f be optically 
thick to neutrino emission, the region of f < r < f be optically thin, and the electron/positron 
pair capture reactions be the dominant cooling mechanism. Thus equation (fTTI) becomes 

f iQjjgSw + f 9 x 10> U T« H2^r = ™ { J_ _ k U _ H (M " 2 1 \ . 
7 r . T Jf 4 [3n r I 3 V r / J J 

(47) 

where we take e ~ 1. Using the self-similar relations and performing some derivations, we 
find 



18.7(1 +Y e f 5 ( m-^m^f-Wj^'^ 



2-7 



r 6 - r *,6 



7-1 



27-1 / 2-37 2-37 
+2.77 ( ^ ) mrridfrj- 1 ( f^ 1 - f^ 1 ) = 0.924mm o 



1+1/(7-1) _ ^,-1+1/(7-1) 
1/2 



.(4f 



The solution shown by equation (f48|) gives f. We take f v ~ 1 and a = 0.1. Moreover, we 
define a new parameter k = f 6 /f 6 , and assume several sets of parameters to give the solution 
of equation ( |48|) . 

From Table 6, we can see that the size of the inner disk increases with the accretion 
rate. In addition, an increase of m or k, or an decrease of Y e also makes the inner-disk size 
larger. We will compare the analytical results from Tables 3 to 6 with numerical results in 
53.4 in more details. 
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3.3. Discussion about the stellar surface boundary 

Now we want to discuss the physical condition near the neutron star surface. We know 
that if the rotational velocity of the neutron star surface is different from that of the inner 
boundary of the disk, then the disk can act a torque on the star at the stellar radius. Here 
we take the rotational velocity of the inner disk Q ~ approximately as mentioned in §2.1 
and §2.3. Then we always have the stellar surface angular velocity to be slower than that 
of the inner disk Qk (i.e., fi* < Qk)- As a result, the kinetic energy of the accreted matter 
is released when the angular velocity of the matter decreases to the angular velocity of the 
neutron star surface. From the Newtonian dynamics, we obtain a differential equation, 

G,r, = <H£L, (49) 
at 

where G*r* = Mr^(Qx — &>*) is the torque acting on the star surface from the disk, I is the 
moment of inertia of the star, I = £,Mr^, with £ being a coefficient. The above equation can 
be further written as 

[n K - + £)]dM = f Mdfl,. (50) 

From equation (|50|) we can see that if Q* t0 < / (1+0 initially, we always have dQ*/dM > 0, 
i.e., the neutron star is spun up as accretion proceeds. On the other hand, if f2* j0 > l+£) 
initially, then the star is always spun down by the disk. The limit value of is 0^/(1 + £) 
in both cases. The solution of equation ([501) is 



M 



M + AM 



£+1 

€ 



n Qk 

1 '*,C 



£ + 1\ AM 



f J M 



1 4*, I 



1+^ 



(51) 

where Mo is the initial mass of the neutron star, and AM is the mass of the accreted 
matter from the disk by the star. The change of f2* depends on the ratio AM/Mq. For 
example, if we assume M = 1.4M , AM = O.O1M and set £ = 2/5, then we have 
(th ~ Q *"> = °- 98 (lTt - Or if we take AM = O.1M , £ = 2/5, and the central stellar 

mass is unchanged, then we obtain (^ — f2*) = 0.79(j^ — fi*,o)- 

Here we consider that the energy released from the surface boundary is also taken 
away by neutrino emission, and assume that neutrinos emitted around the stellar sur- 
face are opaque. We can estimate the temperature of the neutron star surface through 
^aBT 4 Anr^H ~ G ^ r M (2 — e), where H is the half thickness of the inner boundary of the 
disk, and the parameter e has the same meaning as that in §2.3. Then we can estimate 
the surface temperature as T n ~ 0.415(2 — s) 1 l i m 1 l i m l J 4: H^^r^J 2 . This estimation of the 
temperature is only valid if the inner disk is optically thin for neutrinos and only the surface 
boundary is optically thick. If the inner disk, due to a large accretion rate, becomes optically 
thick to neutrino emission, the surface temperature should be higher. 
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3.4. Comparison with numerical results 

In order to give an analytical solution of the accretion disk around a neutron star 
in the simple model, we can choose the dominant terms in equations ffl8|) and ffl~9l) . We 
now consider one type of pressure to be dominated, and assume extreme cases from ADAF 
to NDAF. For example, we consider the neutrino-dominated region with /„ ~ 1 and the 
advection-dominated region with /„ = 0. In this subsection, we solve the hyperaccretion- 
disk structure numerically based on the simple model. In order to compare the numerical 
results with what we have obtained analytically, we keep on with using the equations and 
definitions of all the parameters in this section. However, we first need to point out several 
approximations and some differences between numerical and analytical methods based on 
the simple model. 

First of all, we choose the range of r$ from 1 to 15 in our numerical calculations. In 
other words, we take the range of the accretion disk to be from the surface of the neutron 
star to the radius of 150 km as the outer boundary. During the compact-star merger or 
massive-star collapse, the torus around a neutron star has only some part that owns a large 
angular momentum to form a debris disk, so the mass of the disk may be smaller than the 
total mass of the torus. From view of simulations (Lee & Ramirez- Ruiz 2007), if the debris 
disk forms through the merger of two neutron stars, its outer radius can be slightly smaller 
than the value we give above. On the other hand, the disk size may be slightly larger than 
that we assume above if the disk forms during the collapse of a massive star. The changes of 
physical variables of the disk due to a change of the outer radius may be insignificant, and 
we do not discuss the effect of the outer radius in this paper. 

Second, we fix the viscosity parameter a to be 0.1. If a decreases (or increases), the 
variables of the disk increases (or decreases). More information can be seen in analytical 
solutions in §3.1 and §3.2. In numerical calculations, we take a to be a constant. 

Third, we still set Y e as a parameter in numerical calculations in §3.4. In order to show 
results clearly, we consider two conditions: one is an extreme condition with Y e = 1, which 
means that the disk is made mainly of electrons and protons but no neutrons; the other is 
Y" e = 1/9, which means that the number ratio of electrons, protons and neutrons is 1 : 1 : 8. 
An elaborate work should consider the effect of /3-equilibrium, and we will discuss it in detail 
in §4. 

For numerical calculations, we consider all the terms of the pressure and an intermediate 
case between ADAF and NDAF. In analytical calculations we take the adiabatic index 7 of 
the inner disk to be 5/3 if the disk is gas pressure-dominated or 7 = 4/3 if the disk is radiation 
or degeneracy pressure-dominated. In numerical calculations, however, it is convenient to 



-22 - 



introduce an "equivalent" adiabatic index 7 based on the original definition of 7 from the 
first law of thermodynamics, 7 = 1 + P/u, where P is the pressure of the disk at a given 
r, and u is the internal energy density at the same radius. Therefore 7 is a variable as a 
function of radius. We obtain the self-similar structure, 

P {r) _ t r \ - V(7(r)-1) P ( r ) _ / r \ -7(r)/(7(r)-l) 

p(r + dr) ~ \r + dr ,/ ' P(r + dr) ~ \r + dr ) ' ^ ' 

If 7 does not vary significantly in the inner disk, the difference between the approximate 
solution where 7 is a constant and the numerical solution where we introduce an "equivalent" 
7 is not obvious. 

In addition, some region of the accretion disk is optically thick to neutrino emission 
when rhd is sufficiently large. We use the same expressions of neutrino optical depth and 
emission at the beginning of §3. We also require that the neutrino emission luminosity per 
unit area is continuous when the optical depth crosses r = 1. 

We first calculate the value of f, which is the radius of the boundary layer between the 
inner and outer disks. Figure 1 shows f as a function of for different values of m. We can 
see that when the disk flow is an ADAF at a low accretion rate, f decreases monotonously 
as the accretion rate increases until the value of f reaches a minimum. At this minimum, 
the outer disk flow is an NDAF and most of the neutrinos generated from the outer disk 
can escape freely, and the effect of neutrino opacity is not important. If the accretion rate 
is higher and makes the effect of neutrino opacity significant, the value of f increases with 
increasing the accretion rate. From Figure 1, we see that when the accretion rate is either 
low or sufficiently high, f is very large and even reaches the value of the outer radius, which 
means that the inner disk totally covers the outer disk. In this situation, since the outer 
disk is covered, no part of the disk is similar to that of the accretion disk around a black 
hole as we discussed in §3.1, and thus we say that the entire disk becomes a self-similar 
structure and that the physical variables of the entire disk are adjusted to build an energy 
balance between heating and neutrino cooling. For this situation, we do not want to discuss 
in details any more. We focus on the accretion rate which allows the two-steady parts of 
disks to exist. 

In Figure 2, we choose several special conditions to plot the density, temperature and 
pressure of the whole disk as functions of radius r. Particularly, we take p, P and T to 
be continuous at the boundary between the inner and outer disks. If the parameter m d 
is larger, or the disk contains more neutrons (i.e., Y e becomes smaller), then the density, 
temperature and pressure of the disk are larger, and the change of the density and pressure 
is more dramatic than that of the temperature. The change of the temperature cannot be 
very large because it greatly affects the neutrino cooling rate of the disk. 
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If the mass of the central neutron star (m) becomes larger or the electron fraction 
Y e becomes smaller, then the value of f in the monotonously decreasing segment of — f 
becomes smaller, and the value off in the monotonously increasing segment of rhd—r is larger. 
And the minimum value of f is almost independent of Y e and m. All of these conclusions are 
consistent with the analytical solutions, except for the case of the advection-dominated outer 
disk with the radiation pressure to be dominated. In §3.2.2, we found that the size of the 
inner disk increases with increasing m for our analytical solutions. However, by calculating 
the "equivalent" adiabatic index 7, we find that 7 decreases slightly with increasing m. This 
makes the value of r§ decrease, which is also consistent with the analytical results (see Table 
3). Figure 3 shows the "equivalent" adiabatic index 7 as a function of radius of the entire 
disk for several different sets of parameters. In the case where the accretion rate is low, the 
radiation pressure is important. As the accretion rate increases, the gas pressure becomes 
more dominant and the value of 7 is larger. On the other hand, the ratio of the degeneracy 
pressure to the total pressure is larger in the case of a higher accretion rate and Y e ~ 1. 
However, the gas pressure is always dominant for the accretion rate chosen here. We see 
from Figure 3 that the change of 7 in the inner disk is insignificant. 

Figure 4 shows the ratio of the radial velocity v r and local speed of sound c s as a function 
of r. The ratio is always much smaller than unity, which means that the accretion flow is 
always subsonic and no stalled shock exists in the disk. In many cases the peak of this ratio 
is just at the boundary between the inner and outer disks. The reason can be found in §3.2.1, 
where we gave the analytical expression of the ratio, f = 1 — (7%/r) 1 / 2 is the major factor 
that affects the value of v r /c s of the outer disk. However, in the inner disk, v r /c s always 
decreases since the isothermal sound speed can be greater at smaller radii (i.e., c s oc r -1 / 2 ), 
and the radial velocity of the accreting gas, which satisfies the self-similar solution (TT2"]) . 
cannot change dramatically for the disk matter to strike the neutron star surface. 

Figure 5 shows the total neutrino emission luminosity of the entire disk around a neutron 
st dbT CIS db function of accretion rate for parameters M and Y e (where we do not consider 
neutrino emission from the stellar surface discussed in §3.3), and we compare it with the 
neutrino luminosity from a black-hole disk. Also, we calculate the total neutrino luminosity 
from the inner and outer disks. Here we roughly take the mass of the black hole to be the 
same as that of the neutron star, and the innermost stable circular orbit of the disk has a 
radius which is equal to the radius of the neutron star. We approximately use the Newtonian 
dynamics for simplicity. In Fig. 5, we find that the difference in neutrino luminosity between 
the neutron-star and black-hole cases is a strong function of the accretion rate. When the 
accretion rate is low (m^ < 10), the total neutrino luminosity of the black- hole disk L^bh 
is much smaller than that of the neutron-star disk L^ns, but L^bh and Z^ns are similar 
for a moderate accretion rate (m^ from 10 to 100). Actually, this result is consistent with 



-24- 



the general scenario introduced in §2.1 and the basic result shown in Fig. 1: for a low 
accretion rate, the black-hole disk is mainly advection-dominated with most of the viscous 
dissipation-driven energy to be advected into the event horizon of the black hole, and we have 
L^bh *C GMM /(4r). On the other hand, a large size of the inner disk of the neutron-star 
disk for a low accretion rate makes the neutrino emission efficiency be much higher than its 
black- hole counterpart. However, for a moderate accretion rate, the black- hole disk is similar 
to the neutron-star disk, which owns a quite small inner disk, and we have L^bh ~ -^,ns- 
Moreover, neutrino opacity leads the value of L^bh to be less again compared with L^ns f° r 
a high accretion rate, as this opacity decreases the neutrino emission efficiency in the black- 
hole disk but increases the size of the neutron-star disk again to balance the heat energy 
release. 



4. An Elaborate Model of the Disk 



In the last section we first studied the disk structure analytically. To do this, we used 
several approximations. First of all, we took the pressure as a summation of several extreme 
contributions such as the gas pressure of nucleons and electrons, and the radiation pressure 
of a plasma of photons and e + e~ pairs. However, electrons may actually be degenerate or 
partially degenerate, and the neutrino pressure should also be added to the total pressure. 
Following Kohri et al. (2005), a lot of works about hyperaccretion disks used the Fermi- 
Dirac distribution to calculate the pressure of electrons and even the pressure of nucleons. 
Second, neutrino cooling we used in the last section is simplified, following Popham et al. 
(1999) and Narayan et al. (2001) and neglecting the effect of electron degeneracy and 
the effect of different types of neutrinos and their different optical depth. In fact, these 
effects may be significant in some cases. Third, we took the electron-nucleon-radio Y e as 
a constant parameter in our analytic model in §3. Realistically, Y e should be calculated 
based on /^-equilibrium and neutronization in hyperaccretion disks. In this section, we still 
use the assumption of outer and inner disks discussed in §2, but consider a state-of-the-art 
model with lots of elaborate (more physical) considerations on the thermodynamics and 
microphysics in the disk, which was recently developed in studying the neutrino-cooled disk 
of a black hole. In addition, we compare results from this elaborate model with those of the 
simple model discussed in §3. 
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4.1. Thermodynamics and microphysics 

The total pressure in the disk can be written as: P = P nuc + P ra d + P e + P v . We 
still consider all the nucleons to be free (X nuc « 1) as mentioned in §2.2, and ignore the 
photodisintegration process. Also, we replace the term of radiation pressure llaT 4 /12 in the 
simple model by aT 4 /3 in this section, because the pressure of e + e~ pairs can be calculated 
in the electron pressure P e with the Fermi-Dirac distribution: 

p± _lmy [<*> x* dx 

where x = p/m e c is the dimensionless momentum of an electron and fi e is the chemical 
potential of the electron gas. P e is the summation of P e - and P e +. In addition, we take the 
neutrino pressure to be 

P v = w,/3, (54) 
where u v is the energy density of neutrinos. 

The "equivalent" adiabatic index can be expressed by 

7 = 1 + (-Pnuc + Pad + P e + P V ) / («nue + «rad + U e + U„). (55) 

with the inner energy density to be 

^gas ~^Pgasi (^^) 



Mrad = 3P rad , (57) 
7T 2 h 3 J e (m e c^V^+l T ^)/k B T + l 

We then use equation ( 1521) to obtain the self-similar inner disk. 



In addition, we add the equation of charge neutrality among protons, electrons and 
positrons to estabilish the relation between p, Y e and /i e . 

n p = = n e - - n e +, (59) 

m B 

where we use the Fermi-Dirac form to calculate n e - and n e +. 

Moreover, in the elaborate model, we adopt the improved formula of the neutrino cooling 
rate Q~ , the inner energy density of neutrinos u v , as well as the absorption and scattering 
optical depth for three types neutrinos r a)1/i (e l/i)T ) and r a)1/i ( e)M)T ) following a series of previous 
work (e.g., Popham & Narayan 1995, Di Matteo et al. 2002, Kohri et al. 2005, Gu et al. 
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2006, Janiuk et al. 2007 and Liu et al. 2007). The three types of neutrino cooling rate per 
unit volume are 

— Wefe Qbrcm ~t~ Qplasmonj 

(60) 

Qu M = <lu T = q e -e+^u T Ur + 9brem, (61) 

where the meanings of four terms g e N, Qe-e+^um, <?brem, and <j p iasmon have been shown at the 
beginning of §3. Here g e N and g p iasmon are only related to q Ue . Moreover, g eN is a summation 
of three terms, 

QcN — 4p+e- ->n+u e + Qn+e+^p+u e + <jn-+p+e-+u e ■ (62) 

The formulae of three terms in equation (1621) are the same as Kohri et al. (2005), Janiuk 
et al. (2007) and Liu et al. (2007), who considered the effect of electron degeneracy. In 
addition, we use the same formulae of q e - e +^ u . p ., %rem, and g p i a smon as the early works such 
as Kohri et al. (2002, 2005) and Liu et al. (2007). 

Finally, different from the simple model in §3 in which we took the electron fraction Y e 
as a free parameter, in this section we calculate Y e by considering the /3-equilibrium in the 
disk among electrons and nucleons following Lee et al. (2005) and Liu et al. (2007) 

In (^) = f{T y )^R + [1 - f{r v )p^, (63) 
\n p J k B T k B T 

with the weight factor /(r jy )=exp(— r Ue ) and Q = {m n — m p )c 2 . Equation (163]) is a combined 
form to allow the transition from the neutrino-transparent limit case to the neutrino-opaque 
limit case of the /3-equilibrium. 



4.2. Numerical results in the elaborate model 

Using equations (jSJ), ©, ffTOl and ffTTT) in §2.2 and the improved treatment in §4.1, we 
can solve the structure of the outer disk. Then using equations f|T2l and (IT7|) in §2.3, we 
determine the size and the structure of the inner disk, and calculate the neutrino luminosity 
of the entire disk. 

Figure 6a shows the size of the inner disk in the elaborate model. We still choose the 
mass of the central neutron star to be M = 1.4M and M = 2.OM . From Figure 6a we can 
see that the size of the inner disk r decreases with increasing the accretion rate, and reaches 
a minimum at M ~ 0.1M o s _1 . Then the value of f increases with increasing the accretion 
rate. This result is well consistent with what we have found in the simple model in §3. The 
physical reason for this result has been discussed in §3.3.2 and §3.4. Figure 6b shows the 
solution of the inner disk size both in the simple model and the elaborate model. We fix the 
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central neutron star M = 1.4M . In the simple model of §3, we take Y e as a free parameter 
and plot two rhd — f 6 lines with Y e = 1 and Y e = 1/9, while in the elaborate model we only 
plot one line since Y e can be directly determined through /^-equilibrium in the disk. From 
Figure 6b we conclude that the solutions of the two models are basically consistent with each 
other. 

In Figure 7, we shows the "equivalent" adiabatic index 7 and the electron fraction Y e 
in the entire disk for three values of the accretion rate M = O.OlM s _1 , O.lM s _1 and 
I.OMqS -1 , and for two values of the mass of the central neutron star M = 1.4M and 
2.OM . The "equivalent" adiabatic index 7 increases with increasing the accretion rate in 
most region of the disks, since gas will take over electrons and radiation to be the dominant 
pressure when the accretion rate is high enough. In addition, 7 decreases as the radius 
decreases in the inner disk. These are consistent with the results of the simple model (see 
Fig. 3). From Figure 7b, we can see that Y e ~ 1 when the accretion rate is low, and Y e <C 1 
when the accretion rate becomes sufficiently high. This result is consistent with Kohri et al. 
(2005, their Fig. 6b). Chen & Beloborodov (2007) and Liu et al. (2007) showed the electron 
fraction Y e < 0.5 in the disk, since they supposed that initial neutrons and protons come 
from photodisintegration of a-particles at some large radius far from a central black hole. 
However, since the hyperaccretion disk around a neutron star we discuss has a size smaller 
than that of a black-hole disk, we consider the mass fraction of free nucleons X nuc = 1 in 
the entire disk. Therefore, the fraction of protons can be slightly higher and it is possible 
that the protons are richer than neutrons in the disk or Y e > 0.5 if the accretion rate is low 
enough. 

Figure 8 and Figure 9a show the density, temperature, pressure and the neutrino lumi- 
nosity per unit area in the entire disk with three different accretion rates. We fix M — 1.4M , 
and also plot two other curves of solutions in the simple model with Y e — 1 and Y e = 1/9. 
The density p and pressure P in the elaborate model are smaller than those in the simple 
model when the accretion rate is low (rhd = 1.0), but are similar to the solution of the simple 
model with Y e = 1/9 in the high accretion rate (rhd = 100). In addition, p and P in the 
elaborate model change from one solution (Y e = 1) to another solution (Y e = 1/9) in the 
simple model for an intermediate accretion rate (e.g., rhd ~ 10), since Y e ~ 0.5 in the outer 
edge of the disk and Y e <C 1 in the inner disk. The distribution of the neutrino cooling rate 
Q~ (luminosity per unit area) in the elaborate model is almost the same as that in the simple 
model with Y e — 1 and low accretion rate or Y e = 1/9 and a high accretion rate. However, 
the value of Q~ is still different in these two models for the region that is optically thick to 
neutrino emission in the disk. 

We also plot the the total neutrino emission luminosity of the entire disk, the outer and 
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inner disks as functions of accretion rate with the central neutron star mass of 1.4M Q , and 
compare the total neutrino luminosity with that of the black-hole disk (Figure 9b). The 
results are similar to what we have found in the simple model (Figure 5). 

5. Conclusions and Discussions 

In this paper we have studied the structure, energy advection and conversion, and 
neutrino emission of a hyperaccretion disk around a neutron star. We considered a quasi- 
steady disk model without any outflow. Similar to the disk around a black hole, the neutron 
star disk with a huge mass accretion rate is extremely hot and dense, opaque to photons, 
and thus is only cooled via neutrino emission, or even optically thick to neutrino emission 
in some region of the disk if the accretion rate is sufficiently high. However, a significant 
difference between black hole and neutron star disks is that the heat energy of the disk can 
be advected into the event horizon if the central object is a black hole, but if the central 
object is a neutron star, the heat energy should be eventually released from a region of the 
disk near the stellar surface. As a result, the neutrino luminosity of the neutron star disk 
should be much larger than that in black hole accretion. We approximately took the disk 
as a Keplerian disk. According to the Virial theorem, one half of the gravitational energy 
in such a disk is used to heat the disk and the other half to increase the rotational kinetic 
energy of the disk. We assumed that most of the heat energy generated from the disk is still 
cooled from the disk via neutrino emission and the rotational energy is used to spin up the 
neutron star or is released on the stellar surface via neutrino emission. 

In a certain range of hypercritical accretion rates, depending on the mechanisms of 
energy heating and cooling in the disk, we considered a two-region, steady-state disk model. 
The outer disk is similar to the outer region of a black hole disk. We used the standard 
viscosity assumption, Newtonian dynamics and vertically integrated method to study the 
structure of the outer disk. Since the radial velocity of the disk flow is always subsonic, 
no stalled shock exists in the disk and thus we considered that physical variables in the 
disk change continuously when crossing the boundary layer between the inner and outer 
disks. The inner disk, which expands until a heating-cooling balance is built, could satisfy a 
self-similar structure as shown by equation (TT2l) . 

In this paper we first studied the disk structure analytically. To do this, we adopted a 
simple disk model based on the analytical method. We took the pressure as a summation 
of several extreme contributions and simple formulae of neutrino cooling. And we took the 
electron fraction Y e as a parameter in the simple model. We used an analytical method to 
find the dominant-pressure distribution (Table 2) and the radial distributions of the density, 
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temperature and pressure (solutions [20], ETJ [261 EOl E2]) in the outer disk. Then we used the 
equation of energy balance between heating and neutrino cooling to calculate the size of the 
inner disk in four different cases: whether the advection-dominated outer disk is radiation 
or gas pressure-dominated, and whether the neutrino-cooled outer disk is optically thin or 
thick to neutrino emission (Table 3 to 6). Subsequently, we numerically calculated the size 
of the inner disk, the structure, and energy conversion and emission of the entire disk in 
the simple model (Fig. 2 to Fig. 5) and compared the numerical results with the analytical 
results. The numerical results are consistent with the analytical ones from the simple model. 

When the accretion rate is sufficiently low, most of the disk is advection-dominated, 
the energy is advected inward to heat the inner disk, and eventually released via neutrino 
emission in the inner disk. In this case, the inner disk is very large, and quite different 
from a black hole disk, which advects most of the energy inward into the event horizon. If 
the accretion rate is higher, then physical variables such as the density, temperature and 
pressure become larger, the disk flow becomes NDAF, the advected energy becomes smaller, 
and heating of the inner disk becomes less significant. As a result, the size of the inner disk 
is much smaller, and the difference between the entire disk and the black hole disk becomes 
less significant. Furthermore, if the accretion rate is large enough to make neutrino emission 
optically thick, then the effect of neutrino opacity becomes important so that the efficiency 
of neutrino emission from most of the disk decreases and the size of the inner disk again 
increases until the entire disk becomes self-similar. Besides, a different mass of the central 
star or a different electron-nucleon ratio also makes physical variables and properties of the 
disk different. However, the accretion rate plays a more significant role in the disk structure 
and energy conversion, as it varies much wider than the other parameters. 

The simple model is based on the early works such as Popham et al. (1999) and Narayan 
et al. (2001). We found that the simple model in fact gives us a clear physical picture of 
the hyperaccretion disk around a neutron star, even if we used some simplified formulae 
in thermodynamics and microphysics in the disk. In §4 we considered an elaborate model, 
in which we calculated the pressure of electrons and positrons by using the Fermi-Dirac 
distribution and replaced the factor 11/12 by 1/3 in the radiation-pressure equation. We 
adopted more advanced expressions of the neutrino cooling rates, including the effect of all 
three types of neutrinos and the electron degeneracy. Moreover, we considered /3-equilibrium 
in the disk to calculate the electron fraction Y e . Then, in the elaborate model, we also 
calculated the the size of the inner disk (Fig. 6), the radial distributions of the "equivalent" 
adiabatic index 7, the electron fraction (Fig. 7), the density, temperature and pressure in 
the disk (Fig. 8), the neutrino cooling rate distribution of the disk (Fig. 9a), the neutrino 
emission luminosity from the inner and outer disks, and the total neutrino luminosity of a 
neutron-star disk compared with that of a black-hole disk (Fig. 9b). 
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The electron fraction Y e was also determined in the elaborate model. We found that 
Y e drops with increasing the accretion rate in the outer disk. Y e can be greater than 0.5 
at a large radius if the accretion rate is sufficiently low, and Y e <C 1 in the disk when the 
accretion rate is high enough. If we put these results of Y e in the elaborate model into the 
simple model correctly (i.e. Y e ~ 1 for low accretion rate and F e < 1 for high accretion 
rate), we find that they are basically consistent with each other (see Fig. 6, Fig. 7, Fig. 8, 
Fig. 9), and also consistent with most of the early works (see the discussion in the end of 
§3.1). 

A main difference in the structure between the simple and elaborate models is caused 
by different expressions of pressure adopted in the two models. In order to see it clearly, we 
introduce the third model here. We still keep P ra d = llaT 4 /12 as the simple model but just 
change the relativistic degeneracy pressure term of electrons (the second term in equation 
[T9l) to the Fermi-Dirac distribution (formula [531 for P e -), and use all the other formulae 
in §4. In other words, the third model is introduced by only changing one pressure term 
in the elaborate model. We can find that the results from the third model are even more 
consistent with those of the simple model than the elaborate model in §4. Take Figure 10 as 
an example. We compare the solution of the inner disk of the third model with that of the 
simple model. From Fig. 10 we can see that if the accretion rate is low and Y e ~ 1, the thick 
solid line is much closer to the thin dashed line, which results from the simple model with 
Y e — 1; and if the accretion rate is high enough and Y e <C 1, the solid line is much closer to 
the thick dotted line, which is the result from the simple model with Y e = 1/9. Compared 
with Fig. 8a, this result is even more consistent with the simple model. The values of density 
p and pressure P in Fig. 10, for a low accretion rate, are smaller than those of the simple 
model, which is also due to different expressions of the radiation pressure used in these two 
models in §3 and §4. Therefore, we conclude that the main difference of the results between 
the simple model in §3 and the elaborate model in §4 come from different expressions of 
the pressure in disks. However, we believe that the pressure formulae given in the elaborate 
model are more realistic since llaT 4 /12 is only an approximated formula for the pressure of 
a plasma of photons and e + e~ pairs. On the other hand, as what has been pointed out by 
Lee et al. (2005) and Liu et al. (2007), formulae P ra d = aT 4 /3 and (1531) in §4.1.1 are better 
and can automatically take relativistic e + e~ pairs into account in the expression of P e . 

The different expression of the neutrino cooling rate Q~ makes the neutrino luminosity 
distribution different in the region where is optically thick to neutrino emission. A more ad- 
vanced expression of neutrino cooling rate Q~ gives better results of the neutrino luminosity 
per unit area than that given by the rough expression |<tbT 4 /t in §3. 

In this paper we studied the disk without any outflow, which may exist if the disk flow 
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is an ADAF. However, it is still unclear whether an outflow or neutrino cooling plays a more 
important role since the size of the disk is quite small. The other case that we ignored is 
that, if the radius of the central neutron star is smaller than that of the innermost stable 
circular orbit of the accretion disk, the accreting gas eventually falls onto the neutron star 
freely. In this case, a shock could form in the region between the innermost stable circular 
orbit and neutron star surface (Medvedev 2004). This effect can be studied if other effects 
such as the equation of state of a differentially-rotating neutron star and its mass-radius 
relation are together involved. 

Neutrinos from a hyperaccretion disk around a neutron star will be possibly annihilated 
to electron/positron pairs, which could further produce a jet. It would be expected that such 
a jet is more energetic than that from the neutrino-cooled disk of a black hole with same mass 
and accretion rate as those of the neutron star (Zhang & Dai 2008, in preparation). This 
could be helpful to draw the conclusion that some GRBs originate from neutrino annihilation 
rather than magnetic effects such as the Blandford-Znajek effect. 

We considered a central neutron star with surface magnetic field weaker than B s cr ~ 
10 15 — 10 16 G for typical hyperaccretion rates in this paper. For magnetars (i.e., neutron 
stars with ultra-strongly magnetic fields of ~ -B S;Cr ), however, the magnetic fields could 
play a significant role in the global structure of hyperaccretion disks as well as underlying 
microphysical processes, e.g., the quantum effect (Landau levels) on the electron distribution 
and magnetic pressure in the disks could become important. Thus, the effects of an ultra- 
strongly magnetic field on hyperaccretion disks around neutron stars are an interesting topic, 
which deserves a detailed study. 
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notation 



definition 



j/Eq. 



m 
m d 

Ye 

r* 

r Q ut 

£ 

r 

f 

k = f/r 
/ = 1- 

7=i- 



mass of the central neutron star, m = M/1AM & §3.1, (!20|) 

mass accretion rate, m,i = M/O.O1M s _1 §3.1, (J201) 

ratio of the electron to nucleon number density in the disk §2.2, 

radius of the neutron star §2.3, (TT3"]) 

outer radius of the disk §2.3, (fT3"]l 

angular velocity of the stellar surface §3.3, (}4"9~]) 

efficiency of energy release in the inner disk §2.3, (|17p 

radius between the inner and outer disks §2.3, ( fl2l) 

radius between the neutrino optically-thick & -thin regions §3.2, (|47|) 
parameter to measure the neutrino optically-thick region 

useful factor as a function of r §2.2, ([7]) 

value of / at radius f §3.2, (13"5|) 

average efficiency of neutrino cooling in the outer disk §2.3, fll3p 



Table 1: Notation and definition of some quantities in this paper. 
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Table 2: Range of the accretion rate in different regions, (a) rh d < 0.504m 1 / 2 a^ 1 3 r , g //2 i /' 1 / 6 ; 
(b) m d < 644(1 + Y e y 4 m 7 / 2 a_ ir - 3/2 f 7 / 2 ; (c) m d > 644(1 + Y e )- 4 m 7 / 2 a^r 6 3/2 fl 2 - 
(d) m d < 4.81 x 10- 3 (1 + y e ) 13 / 5 m- 29 / 2 VY 10 ^f /2 °r 2 ; (e) m rf > 4.81 x 10~ 3 (1 + 

y e) 13/5 m -29/20 a 21/10 r 47/20 / _ 2 . (f) ^ < Q 14g y-4 (1 + y j27/5 m -ll/10 a 19/10 7 .33/20 / _ 1 . (g) 

m, > 0.148y~ 4 (l + Y e ) 2r / 5 m- n / 10 a^{ 10 4 3/20 f- 1 . 



m d 


0.01 


0.05 


0.1 


0.5 


Case 1 


6.38 


3.53 


2.76 


1.65 


Case 2 


13.87 


5.66 


3.92 


1.84 


Case 3 


6.83 


3.76 


2.94 


1.74 



Table 3: Equation ( [4T|) gives r$ in several cases. Case 1: m = 1,7 = 4/3; case 2: m = 1,7 = 
1.4; case 3: m = 2.0/1.4, 7 = 4/3. 
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m d 0.5 1 3 5 10 

Case 1 2.17 1.91 1.62 1.52 1.42 

Case 2 2.05 1.83 1.59 1.50 1.40 

Case 3 1.68 1.54 1.38 1.33 1.27 

Case 4 1.97 1.76 1.52 1.44 1.36 



Table 4: Equation (145]) gives r$ in several cases. Case 1: Y e = 1,7 = 5/3, m = 1; case 2: 
Y e = 1,7 = 3/2, m = 1; case 3: Y e = 1/9,7 = 5/3, m = 1; case 4: Y e = 1,7 = 5/3, m = 
2.0/1.4. 



M7 


5/3 


3/2 


4/3 


0.9 


1.31 


1.28 


1.26 


0.7 


1.46 


1.43 


1.39 


0.5 


1.56 


1.53 


1.47 



Table 5: Equation ( l46l) gives r§. We take f v and 7 as parameters to solve this equation. 



m d 80 100 120 150 

Case 1 4.28 5.13 5.94 7TlT 

Case 2 6.26 7.48 8.65 10.33 

Case 3 5.09 6.09 7.05 8.43 

Case 4 4.10 4.99 5.86 7.11 

Case 5 4.87 5.90 6.89 8.32 



Table 6: Equation fj4"8l) gives r 6 in several different cases. Case 1: Y e — 1, m — 1, k — 1, 7 = 
5/3; case 2: Y e = 1/9, m = l,k = 1,7 = 5/3; case 3: F e = l,m= 2.0/1.4, fc = 1,7 = 5/3, 
Case 4: F e = 1, m = 1, = 0.7, 7 = 5/3; case 5: Y e — 1, m = 1, k — 1, 7 = 3/2. 
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Fig. 1. — The radius r$ of the boundary layer between the inner and outer disks as a function 
of accretion rate rhd in the simple model for several sets of parameters: (1) M = 1.4M and 
Y e = 1.0 {solid line), (2) M = 1.4M and Y e = 1/9 (dashed line), (3) M = 2.OM and 
Y e = 1.0 (dotted line), and (4) M = 2.OM and Y e = 1/9 (dash-dotted line). 
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Fig. 2. — The density (in units of 10 11 g cm -3 ), pressure (in units of 10 11 ergs cm -3 ) and 
temperature (in units of 10 n K) of the disk as functions of radius r (in units of 10 6 cm) in the 
simple model for several sets of parameters: (1) M = 1.4M , Y e = 1.0, and M = 0.01M o s -1 
(thin solid line), (2) M = 1.4M , Y e = 1.0, and M = O.lM s -1 (thin dashed line), (3) 
M = 1.4M , Y e = 1.0, and M = l.OM s -1 (thin dotted line), (4) M = 2.OM , Y e = 1.0, and 
M = O.OlMoS" 1 (thick solid line), (5) M = 2.OM , Y e = 1.0, and M = O.IMqS' 1 (thick 
dashed line), and (6) M = 2.OM , Y e = 1.0, and M = l.OM s -1 (thick dotted line). 
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Fig. 3. — The "equivalent" adiabatic index 7 of the disk as a function of radius r in the 
simple model. The meanings of different lines are the same as those in Fig. 2. 
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Fig. 4. — The ratio of the radial velocity v r and local speed of sound c s as a function of r. 
(a) Left panel: M = 1.4M , Y e = 1.0; (b) Right panel: M = 1.4M , Y e = 1/9. The accretion 
rate: m d = 0.01M o s _1 (solid line), m d = 0.1M o s _1 (dashed line), and fn d = l.OM s _1 
(dotted line) 
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Fig. 5. — Neutrino luminosity from the disk around a neutron star in the simple model. 
L v is in units of 10 51 ergss _1 . (a) Left panel: M = 1.4M and Y e = 1.0. (b) Right panel: 
M = 1.4M and Y e = 1/9. The solid line corresponds to the maximum energy release rate 
of the disk around a neutron star, the dashed line to the neutrino luminosity from the inner 
disk, the dotted line to the neutrino luminosity from the outer disk, and the the dash-dotted 
line to the neutrino luminosity from a black hole disk. 
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Fig. 6. — (a) Left panel: the radius f 6 between the inner and the outer disks in the elaborate 
model with M = 1.4M (solid line), and M = 2.OM (dashed line), (b) Right panel: 
comparison of f 6 in different models with M = 1.4M : (1) the elaborate model (solid line), 
(2) the simple model with Y e = 1 (dashed line), and (3) the simple model with Y e = 1/9 
(dotted line). 




Fig. 7. — (a) Left panel: the "equivalent" adiabatic index 7 in the elaborate model, (b) Right 
panel: the electron fraction Y e in the elaborate model as a function of radius. The profiles are 
shown for three values of the accretion rate and two values of the central neutron star mass: 
(1) M = 1.4M and M = O.O1M s" 1 (thin solid line), (2) M = 1.4M and M = O.1M s" 1 
(thin dashed line), (3) M = 1.4M and M = l.OMgS" 1 (thin dotted line), (4) M = 2.OM 
and M = O.OIMqS^ 1 (thick solid line), (5) M = 2.OM and M = O.lM s _1 (thick dashed 
line), and (6) M = 2.OM and M = l.OM s' 1 (thick dotted line). 
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Fig. 8. — The density (in units of 10 11 g cm" 3 ), pressure (in units of 10 11 ergs cm -3 ) and 
temperature (in units of 10 n K) of the entire disk in two models for M = 1.4M . The profiles 
include three groups of lines and each group also includes three lines, which are shown for 
three values of the accretion rate M = O.OlM s _1 , O.lM s _1 and l.OM s _1 from bottom 
to top in these figures. The solution in the simple model with Y e = 1 is shown by the thin 
dashed line, and Y e — 1/9 by the thin dotted line. The solution in the elaborate model is 
shown by the thick solid line. 
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Fig. 9. — (a) Left panel: the neutrino luminosity per unit area (in units of 10 39 ergs cm -2 s _1 ) 
in both the simple model and the elaborate model. The meanings of different lines are the 
same as those in Fig. 8. (b) Right panel: the neutrino luminosity from the disk in the 
elaborate model. The meanings of different lines are the same as those in Fig .5. 
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Fig. 10. — Comparison of f 6 in three models with M = 1.4M Q : the simple model with Y e — 1 
(dashed line), the simple model with Y e = 1/9 (dotted line), and the third model discussed 
in §5 (thick solid line). The solution of f 6 in the third model is even more consistent with 
the simple model than that of the elaborate model discussed in §4. 



